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Abstract 

Gallavotti proposed an equivalence principle in hydrodynamics, which states that forced-damped fluids can 
be equally well represented by means of the Navier-Stokes equations and by means of time reversible dynamical 
systems called GNS. In the GNS systems, the usual viscosity is replaced by a state-dependent dissipation term 
which fixes one global quantity. The principle states that the mean values of properly chosen observables are 
the same for both representations of the fluid. In the same paper, the chaotic hypothesis of Gallavotti and 
Cohen is applied to hydrodynamics, leading to the conjecture that entropy fluctuations in the GNS system 
verify a relation first observed in nonequilibrium molecular dynamics. We tested these ideas in the case of two- 
dimensional fluids. We examined the fluctuations of global quadratic quantities in the statistically stationary 
state of a) the Navier-Stokes equations; b) the GNS equations. Our results are consistent with the validity of 
the fluctuation relation, and of the equivalence principle, indicating possible extensions thereof. Moreover, in 
these results the difference between the Gallavotti-Cohen fluctuation theorem and the Evans-Searles identity is 
evident. 

PACS1998 numbers: 47.27.Gs, 05.40. +j, 05.70.Ln 
AMS1991 numbers: 82C05, 76F20 
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> '■ 1 Introduction 

The fundamental laws governing the behaviour of fluids are known, and universally believed to be correctly repre- 
C$ ' sented by the Navier-Stokes (NS) equations. However, a clear cut connection between microscopic and macroscopic 
scales, which definitely justifies such a belief, has still to be made, particularly for systems subjected to nonconser- 
vative external fields Q. If one considers this kind of problems as merely technical, and accepts the NS equations 
as a valid tool for the description of a fluid's dynamics, serious difficulties are met nevertheless in the study of 
the mathematical properties of such equations. For instance, the problem of global existence and regularity of 
solutions of the NS equations is far from being solved in as general terms as desired. This fact is not of exclu- 
sively mathematical interest, because it casts doubts on the effectiveness of the approximation methods devised 
to extract information from the NS equations. Therefore, what is reasonable and presently feasible is to study 
simplified problems, and to construct theories connecting such problems with the full NS equations or with the 
real dynamics of the fluids at hand. In Ref.Q, this is stated as: "It has been recognized that a realistic goal for a 
statistical theory of turbulence is to determine the equations governing the dynamics of some reduced set of modes 
which allow calculation of fundamental quantities of the flow". 

The situation is similar, to some extent, to that of molecular dynamics, in nonequilibrium statistical mechanics. 
There, some progress took place when infinite reservoirs, or driving boundary conditions were replaced by artificial 
constraints imposed on the bulk dynamics of TV-particle systems, with N < oo (see e.g. |], |^, ^|). This way direct 
numerical simulations of the particle models become feasible, and dynamical systems theory leads to theoretical 
predictions Q which can be tested in numerical simulations or real experiments. The trade-off of this approach is 
that the relevant dynamical equations do not seem to be fully justified on physical grounds: no known fundamental 
force acts on the particles of the system in such a way as to implement the desired constraints. Nonetheless, the 
results obtained this way are in excellent agreement with experience [|| , and several arguments have been developed 
to explain why this should be the case. Among such arguments, we deem more convincing those invoking a kind 
of equivalence of ensembles, which is a well known concept in equilibrium statistical mechanics, but rather new 



in noncquilibrium statistical mechanics @, If verified, the equivalence of ensembles guarantees that different 
microscopic dynamics result in the same macroscopic behaviour, thus justifying the use of one kind of dynamics or 
another, depending on which one is more natural to study a given physical problem^] 

Similarly, in Ref.||, Gallavotti argues that different probability distributions characterizing the steady state 
properties of a fluid should yield the same values for (some, at least) macroscopic quantities, if such distributions 
are obtained from different microscopic models of the same fluid. In particular, inspired in part by the mentioned 
developments in statistical mechanics jfj] and by Ref.[||, Gallavotti conjectures one kind of equivalence between the 
NS equations and special time reversal invariant equations which he calls GNS, where the G stands for Gaussian 
constraint. 

The paper (|] introduces a constrained Euler system, which encompasses a portion of the degrees of freedom 
of the NS equations sufficient to obtain some of the fundamental statistical properties of the fluid. The explicit 
form of the constrained equations is derived from the observation that in stationary isotropic turbulence, the mean 
energy in a narrow wave-number shell is nearly constant in time. These equations are similar to those used in 
nonequilibrium molecular dynamics, for driven particle systems subjected to a "gaussian thermostat" || [|. 
Therefore, the properties of such particle systems should be observed to some degree in the constrained Euler 
system of || and, perhaps, in still more general settings. This led Gallavotti § to revive and put under a new light 
Ruelle's principle for hydrodynamics p0[ . This principle had been previously extended to nonequilibrium statistical 
mechanincs in ||, where the first proof of the Gallavotti-Cohen Fluctuation Theorem (GCFT) was obtained, under 
the assumptions of the Chaotic Hypothesis (CH) quoted below in Section 2. 

Tests of the validity of the equivalence conjecture (EC) and of the CH include a Raylcigh-Bcnard convection 
experiment by S. Ciliberto and C. Larochc |jjfl , and numerical simualtions of the GOY shell models by L. Biferale, 
D. Pierotti and A. Vulpiani p^| . The experiment of |Tl[ is consistent with the validity of a fluctuation relation || 
similar to that of the GCFT; while the simualtions of |12[ evidenced a kind of equivalence of different hydrodynamic 
models: the equivalence of energy cascades. A posteriori, the work by She and Jackson j|] can also be taken as a 
verification of the validity of the EC, although apparently not motivated by a general theory. Among other relevant 
works, notable are Refs.gl, H H, although they do not concern hydrodynamic equations. 

In the present paper, we test both the validity of the mentioned fluctuation relation and of the EC for the 
NS and GNS equations. Because lengthy calculations are needed, we consider two-dimensional systems rather 
than three-dimensional ones. Our results confirm, and actually indicate possible extensions of both the EC and 
the fluctuation relation of Ref . [|| , consistently with Gallavotti's predictions for the slope of such (linear) relation. 
Furthemrore, as explained in Section 5, these results provide an example in which the difference between the GCFT 



and an identity previously obtained by Evans and Searles 16 is evident 



2 Reversible damping and the equivalence conjecture 



The CH has been introduced 
and to "gaussian" constraints 



[S| in the study of Af-particle systems subjected to nonconservative external forces 
\2 , ||, D, like dynamical systems of the form 
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defined on a (2dN — l)-dimensional manifold £1 C M 2dN , d being the dimension of the physical space. Here, (q^, p;) 
is the usual notation for position and momentum of particle i, F m ^ represents the action of the other N — 1 particles 
on particle i, Cj is a charge coupling particle i to the external field F ex ^, and ap.i is a dissipation term, which allows 
the system to reach a stationary state. In particular, the choice 
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where K = X)iliP?/2 m i is the kinetic energy of the system, implies that the internal energy of the system 
E = K + $, $ being the internal potential energy, is a constant of the motion. We denote by Vt : ^ — » f2, 
t € M, the flow, so that t \— > Vtj represents a solution of Eqs.Q with initial condition 7. The dynamical system 
(Vt, fi) is but an example from a wide class of systems which share remarkable properties |(|. For example, such 
systems are time reversal invariant, although dissipative. This means that the involution i : {pj}^^ <— > {~Pj}jLi 
anticommutes with the time evolution: iVt — V-t% (reversibility); but phase space volumes contract on average 
(dissipativity), resulting in a multi- fractal stationary state ||l7f . Such systems have been studied in detail, under 
the assumption that the following holds §: 

1 Consider, for instance, the success of the idealized dynamics of lattice gas cellular automata in describing quite complex hydrody- 
namic situations tel. 



Chaotic Hypothesis (CH). A chaotic many-particle system or fluid in a stationary state can be regarded, for 
the purpose of computing macroscopic properties, as a smooth dynamical system with a transitive Axiom- A global 
attractor. In reversible systems it can be regarded, for the same purpose, as a smooth transitive Anosov system. 

This approach led to interesting results for the macroscopic properties of given systems, directly from their mi- 
croscopic dynamics. Among these are proofs of the positivity of transport coefficients jl8[ of the validity of 
Onsager relations^] [^2| and of the fluctuation relation of 0, whose hydrodynamic version Eq.(p6|) is given 
below. 

Following Gallavotti's ideas ||, we now consider the NS equations for a newtonian incompressible fluid: 

ii+ (u • V)u = --\7p + g + i^iAu , V-u = 0. (3) 

P 

Here, u is the velocity field, p is the fluid density, p the pressure, g is a constant forcing term and v\ is the constant 
viscosity. The curl of this equation gives 

ui + (u ■ V)uj = (lj ■ V)u + f + (/[Aw (4) 

in which to = V x u is called vorticity and f = V x g represents the forcing term. If we replace V\ by 

J [w ■ f + u ■ (to ■ V)u] dx 
J (V x u>) 2 dx 



/3 (u,w,f) = =£-l r/ „ V x2 / J , (5) 



we obtain a system whose total enstrophy Q = J ui 2 dx is a constant of motion, and which is time reversal invariant 
in the sense Eqs.(|l|) are. Similarly to Ref.[|), we refer to Eq.(||) with f3q in place of v\, as to the GNS equations, 
i.e. the NS equations with a Gaussian constraint. In two spatial dimensions, the vortex-stretching term {to ■ V)u 
vanishes; hence only a dynamical equation for the third component of the vector uo is necessary. Denoting this 
component also by lo, the vorticity equation reduces to 

lu = -(u- V)uj + f + aAuj. (6) 

where a stands either for v\ or for @q, depending on the case. We impose doubly-periodic boundary conditions 
on Eq.(|^), and rescale the size of the system to 27r, making natural a Fourier expansion for lu. Such expansion, 
truncated for numerical implementation, yields the approximation 

N N 

«(x)w J2 J2 elk ' Xw k, k=(k x ,k y ), NtN. (7) 

k^ = -N k y =-N 

Substituting Eq.(^) in Eq.(|^), we get the equations for the Fourier modes w^. We write such equations in the more 
general form: 

w k -r k + / k -ak 2! w k , k x ,k y = -N,...,N , Z = 0, 1,2, ... (8) 

where each I accounts for the dissipation given by a power of the Laplace operator applied to to. This way, we get 
Ekman damping for I = (a — vq), normal viscosity for I = 1 (a = fx), and a different hyperviscosity vi for any 
I > 1, while only / = 1 was considered in H. The quantity rk is the nonlinear interaction term, which reads 



Equations (||) determine the dynamics of the (2N + l) 2 complex modes w k , in a phase space Q C C < - 2A,+1 ' whose 
dimension is only 2N(N + 1), because of the reality condition lj^ = cj^ k and of the absence of the k = (0, 0) mode. 

Differently from we consider all the quadratic global quantities, such as the energy E — J wA^^dx, the 
enstrophy Q = J w 2 rfx, the palinstrophy P — J (Vw) 2 dx, the "hyperpalinstrophy" H = J (Alj) 2 g?x, etc, instead of 
considering E and Q only. Using the spectral notation, we can write any one of them as 

Q m = ^Tk 2 "V^k, (10) 



2 The microscopic models considered in these studies constitute a wide class with some peculiar feature, such as the unusual dissipation 
terms. Therefore, it is not immediately clear that locutions borrowed from physics —such as "Onsager relations"- can be attributed 
the same phenomenological meaning they_usually have. For instance, a minimal prerequisite for the models to be physically relevant 
is that their number of particles be large [20| , while this is not assumed in Jig - fell] . 



where Q-\ = E, Qo = Q, Qi — P, Q2 — H etc. These quantities satisfy the evolution equations 



where we have defined for convenience of notation 
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and I depends on the kind of viscosity under consideration. The term R m fi vanishes for m = — 1 and m = 0, 
because of the symmetries of the summands: this corresponds to the conservation of energy and enstrophy in 
inviscid two-dimensional flow. Now, if we take 



:/3, 
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(13) 



instead of the constant viscosity coefficient vi in Eq.(J8|), Q m becomes a constant of motion. In this case, we refer 
to Eqs.(||) as to the "cut-off GNS equations". [] An important quantity is a a , the negative of the divergence in 
phase space of the right hand side of Eqs.(||), 
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which may be identified |22j, M with the entropy production rate of the system. Then, for a = vi we get 
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fe x = -AT k v =-N 



while for a = /3; >m (the GNS case) we have 

„GNS _ a 
a l,m - Pi,: 

The sums in Eq.([l5|) evaluate to 
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/ = (Ekman damping), 
I = 1 (normal viscosity), 

I = 2 (hyperviscosity), etc. 



(17) 



Clearly, af s is a constant, while crf^ 3 depends on the vorticity field, hence on time. Nevertheless, /3; jm could fluc- 
tuate around an average value (pi, m ) in such a rapid fashion as to be unresolved on macroscopic scales, making the 
macroscopic behaviour of the GNS model identical to one with given constant viscosity. This is our interpretation 
of the EC of H which for I = 1, m — states that: 

Equivalence Conjecture (EC). The stationary probability distributions of the NS equations and of the GNS 

equations are equivalent in the limit of large Reynolds number, provided the enstrophy Q and the viscosity v are so 
related that o~f s and the average of o~f^ s are equal. 

the large Reynolds number plays the role of the thermodynamic limit, in which limit Pi tm o"f s / vi in 



is expected to dominate over the other terms, making the equality af s 



(of*" 



) imply (A,m) = v u for 



In Ref.| 
Eq.@ 

I = 1, r7i = 0. We further conjecture that the same should hold for any I and m. 

How can these reasonings be justified? In the first place the solutions to the NS equations are expected to take 
a form, in the vi — > limit, which is independent of the viscosity and of the boundary conditions, and therefore 
universal. This is because the flow is expected to develop on characteristic length scales which are much smaller 
than the system size but much larger than the small scale at which the dissipation is effective. This limit is 
then natural in hydrodynamics, as it is desirable for the EC to have a universal rather than particular character. 



3 Equations (g), with a = /3; im , for an arbitrary choice of l,m, cannot be interpreted in terms of Gauss' principle, while that 
interpretation was allowed in the cases of Ref. |i| . 



Moreover, in standard computational practice, the limit vi — » is connected to the N — > oo limit: a numerical 
simulation of turbulence is considered reliable when the dissipative scale 1/kd is well resolved (iV > Lkd, L being a 
length scale of the system). The value kd is estimated by applying Kolmogorov's theory, which assumes a power law 
scaling in a range of k called inertial; this scale becomes smaller as vi — > 0. However, the GNS case is more garbled, 
as it is not known how the scaling is affected by the truncation and by the reversible dissipation. Arguments based 
on power counts shed little light on this, at least for the easier cases in which R m d = 0: for increasing N, <r„ s 
behaves as 

A x rN N 2l+2 _ 1 



k T =-N fe„=-AT 



^ E k a, ~2M7r/ fc 2,+1 rffc~n 2l + 2 . (18) 



In order to derive scaling relations in TV for the quantities involved in the definition of uf^ 8 , we can assume that 
the shell energy E(k), defined as 

£(fc)=£^, (19) 
[|k|]=fe 

where [•] denotes the nearest integer, scales in k as k x . Such a power law provides at least a standard term of 
comparison. The reference theory for 2D turbulent flows, due to Batchelor, Kraichnan and Leith pa], predicts 
A = —3. This slope would appear only in the inertial range, and a steeper decay of E{k) is expected in the 
dissipative range k > kd- It can then be inferred that, for N large, but smaller than Lkd, Qn scales in N as 



N N r N AT2n+A+3 _ 1 

Qn = E E k2 >x| 2 - E k 2n+2 E(k) ~ 2. / k^^dk ~ ^ 1 . (20) 

k u: =-Nk v = -N k Jl 

This implies that Q n converges for large N if A < 2n + 3, while it would diverge in the opposite case. This 
divergence might indeed not be observed, as the limit N — > oo has to be taken before v\ — ► 0, i.e. a dissipative 
range is always included. As for F m in Eq. (|l3|) , we argue that for a given forcing concentrated on a few (or just one) 
k vectors, the term F m can depend on N only if the magnitudes of |o>k| or the phase between uj^ and /t change 
when adding more ultraviolet modes. If we assume that the (large scale) vorticty field is not significantly altered 
by the (much smaller) additional large N modes, then also F m should not change significantly while increasing N. 
Eq.(|l3|) therefore says that 

Pl ' m ~oL~ \N^+^-l\ > m =~ 1 >°- ( 21 ) 

Our results reported in Section 4 show that (3i :in ~ vi, for all N, so we should have A < —2m — 21 — 3 (e.g. A < —3 
for constrained energy and A < —5 for constrained enstrophy and normal viscosity), provided that E{k) ~ k x and 
within the inertial range. With this constraint for A, we have 



Chi+rn _ ( 0(N 21 ) if > A + 2m + 2Z + 3 > -21 



Ql+m \ O(N ) if A + 2m + 21 + 3 < -21 
The power count in eq. ([l6]) then shows that for I > we have 



(22) 



ofNS = ^L 0(7V 2/ +2) ^ (23) 
Vl 

so that for large N (depending on v\) af^ s is proportional to /3z, m , in the leading order. Indeed this is correct for 
any N, if / = 0, and numerical inspection shows that, in all our simulations (made with / = 1 and m = — 1, 0, 1, 2), 
the ratio vii&f 3 ^ 3 ) / {erf 18 (/3;, m )) differs from one by no more than a few parts per thousand, already at truncations 
with N as low as 3. We can therefore fix values of vi and number of modes N, and check the equality (/% >m ) = V\ 
needed for the macroscopic equivalence of the NS and GNS systems, from which the equality af s — (o~f ! ^ s ) 
follows. 

Whether we take the thermodynamic limit or not, there is a formal similarity between the cut-off GNS equations 
and Eqs. (]l]), which may then enjoy similar properties. Of these, the validity of the GCFT of || is perhaps the 
most striking one, and in || Gallavotti has argued that a similar relation (introduced in |13]]) should hold for the 
fluctuations of o~f^ s , defined by Eq.(|T6|) with / = 1 and m = 0. Because of Eq. (p3|) , a similar fluctuation law 
should also hold for the fluctuations of [3i >m (possibly for any I, m), which has a more direct physical interpretation 
than crfm S i sm ce /3; iTn represents the viscosity in the GNS equations. 

To test these ideas in our simulations of a 2D fluid, we need the following definitions. Let 1 1— > Vtcu be the time 
evolution for an initial vorticity field lo, take an overall simulation time T which must be adequately longer than 



the characteristic time of the fluctuations of /3i iTn , and assume that relaxation to a statistically stationary state 
takes place in a time to <C T. Then, introduce the time average 



T 



(Pi,m) = y / Pi, m (Ytu)dt. (24) 







We omit to indicate the dependence of {/3i <m ) on T and on u), assuming that T is large enough, and the dynamics 
ergodic enough, that larger simulation times and different randomly chosen initial fields would produce negligible 
differences in the end result. Let us subdivide the time interval [to, T], in a number of subintervals of length r, and 
consider the quantities 



'','„,(') = — / ft, m (V t w) dt ; i = l,...,^J±. (25) 

/t +(i-l)T T 
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These values, arranged in a histogram, allow us to construct the probability distribution ir[ m of the fluctuations 
of /3/, m integrated over a time r. Then, a test of the validity of the fluctuation relation Eq.(4.1) of |§ amounts in 
our setting to verify that the values 

obtained for different values of r and p (in the support of 7r ; r m ), do converge for growing r to a definite value Co n 
independent of p. The fluctuation relation of || is then a statement on the form of the measured p.d.f. 7r ; r TO for 

large r. We found that C l ^ m (T,p) can be fitted by a p-independent function of r, C^ m (r) say. Therefore, the odd 
part of log(7r ; r m (p)) appears to be linear in p at all r. Then, to assess the validity of the fluctuation relation, we 

verify whether the slope CL' m (t) of the linear function of p 

1 



T~{Pl, m ) 



N « m (p)) - log « W (-P))] = C l f{r)p (27) 



does converge to a given c^ m , as r is increased. This linear form is shared by many kinds of distributions, including 
distributions with exponential tails and Gaussian distributions, giving the erroneous impression that a trivial 
connection exists between the fluctuation relations of || ^ and the central limit theorem (see point 4 of Section 
5). 



3 Description of the numerical method 

We performed a number of numerical experiments, integrating in time several vorticity fields and using a standard 
protocol for the NS equations |Q. We used a two-dimensional 2/3 dealiased pseudospectral code on the periodical 
square [0,27r] x [0,27r]. The resolution of such simulations is traditionally expressed by the number of modes 
M X M over which the convolution term is computed; the number M is often a power of 2 for FFT convenience. 
With the dealiasing procedure, the vorticity field is reconstructed in physical space using only the modes with 
—M/3 < k x ,k y < M/3. The formulas given above have therefore to be read with N = M/3: for example, 
M = 32 implies a system of 221 complex coupled modes. Time advancement in our algorithm is accomplished by 
a fourth order Runge-Kutta integrator. This NS code is easily modified to integrate the GNS equations; even the 



convolution term R m .i in (12) can be computed via a multiplication in physical space, as is done according to the 
pseudospectral method for the nonlinear term rife of eq. (|9j). However, in the GNS case the analytical integration of 
the linear viscous term, often done to allow longer integration steps for the NS equations, is not possible. Moreover, 
the mere substitution of the constant vi by /3z, m , to be computed at each substep of the Runge-Kutta algorithm, 
does not guarantee the exact conservation of the constrained global quantity Q m . This is not a real problem, 
because the relative error e = \Q m (to) — Q m (to + T) \ /Q m {to) made in the conservation of Q m can be controlled 
by an adequate reduction of the integration time steps. In Table |l|, the values of e for each run are explicitly given, 
confirming that our accuracy in the conservation of Q m is satisfactory. 

To study the properties of statistically stationary, nonequilibrium states both for the NS and for the GNS 
systems, we need a criterion to judge if a stationary state has been reached in a simulation. We found that, for the 
NS system, the approach to a steady state can be assessed by looking at the time series of the quantity Q m which 
is going to be constrained in the GNS system. Fluctuations of such quantities in the NS system have amplitudes 
and correlation times which vary from case to case, but they take relatively short times to settle around an average 




Figure 1: Run #9: behaviour of E, Q, P in the NS "buildup" part (first halves of the plots), followed by the 
behaviour of the same quantities and of once E is constrained according to the GNS dynamics (second halves 
of the plots). Data points are subsampled for graphical pourposes. {Pi-i) is close to v\. 



value, if vi is not too small and the forcing is not too large.Q As an empirical criterion, we then assume that the 
NS system has reached a stationary state when Q m fluctuates around a definite value. Figures |l| and ^| show how 
this is typically achieved. 

To perform our NS simulations, we have first to choose the values of the viscosity coefficient vi, of the force /k, 
and of the number of active modes, i.e. of N. The correspondig GNS system is obtained keeping the same /k and 
N, discarding vi, and choosing the global quantity to be fixed which, in our simulations, was E or Q or P or H. 
Once this is done, we let the NS system evolve towards a steady state using a "buildup" approach: starting from 
a small random initialization of the modes, we numerically integrate the field for a sufficiently long time, t 
say. Simulations are considered sufficiently long in this context, when they are long with respect to the correlation 
times of the fluctuating series and to the time required to stabilize the average of the considered Q m . A snapshot 
of the vorticity field at an arbitrary time t > to, is then used as initial datum for the GNS runs. 

In order to perform longer runs, we minimized the number of modes of the vorticity, i.e. the resolution. Some 
numerical experiments led us to conclude that, given a force /k, an adequate resolution should include the spectral 
modes up to at least 2k. Thus, for instance, we performed many of our runs with forcings on k = (±3, ±4) and 
k = (±4, ±3) at a resolution of M = 32, which, taking into account the dealiasing, implies a truncation of the 
Fourier expansions at N = 10. Such resolution allows us to carry on simulations for up to some million timesteps, 
i.e. to total simulation times T of the order of 10 5 time units, whereas fluctuations show characteristic times of 
the order of ten units. Runs at higher resolution converge to a steady state with approximately the same values of 
the chosen Q m , in roughly the same time and, in general, with similar fluctuation amplitudes. The values of the 
Reynolds numbers (computed as are n °t particularly relevant in our context; however, they are observed to 

be of the order of 10 2 up to 10 4 , in the various runs. The complete list of runs performed, inclusive of simulation 
parameters and results, is presented in Table fil. 



4 Tests: equivalence conjecture and fluctuation formula 

Our GNS systems, with our choice of parameters, take a reasonable computation time to reach a stationary state, 
i.e. a state in which the fluctuations of /3; jm occur about its time average (/3z,m). This quantity is found to coincide 
(apart from small errors, and from the cases with insufficient statistics) with the value of the viscosity vi of the NS 
equations, if the quantities E, Q. P and H are fixed to the value measured in the final NS snapshot. This fact is 

4 This and the fact that sufficient numerical accuracy requires very small integration time steps if ui is small and / large, determines 
the range of parameters which we can investigate. 
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8 


256 


0.1 


5 


p 


2 • 10" 6 


10/55 


1 


1.6024 




9 


32 


0.025 


150 


E 


5 • 10" 5 


IOO/55 


1 


0.9497 


0.30 ±0.01 


10 


32 


0.025 


150 


Q 


4 • 10~ 3 


IOO/55 


1 


0.8891 


0.35 ±0.03 


11 


32 


0.025 


125 


p 


1 ■ 10~ 2 


IOO/55 


1 


0.9344 


0.27 ±0.04 


12 


64 


0.025 


25 


E 


9 • 10~ 6 


IOO/55 


1 


1.0189 


0.011 ±0.001 


13 


64 


0.025 


25 


Q 


5 • 10~ 3 


IOO/55 


1 


1.0359 


0.0188 ±0.0006 


14 


64 


0.025 


25 


p 


3 • 10- 2 


IOO/55 


1 


0.9837 


0.075 ±0.01 


15 


32 


0.002 


2 


E 


5 • 10- 8 


1000/55 


1 


0.9439 


0.37 ±0.03 


16 


32 


0.002 


2 


P 


1.5- 10" 5 


IOOO/55 


1 


1.02558 


0.19 ±0.02 


17 


32 


0.025 


150 


E 


8 ■ 10" Y 


IOO/55 


3 


0.9406 


0.23 ±0.05 


18 


32 


0.025 


150 


Q 


2 ■ 10" 4 


IOO/55 


3 


0.7553 


0.36 ±0.08 


19 


32 


0.025 


150 


p 


9 • 10" 6 


IOO/55 


3 


1.6349 


- 0.2 


20 


32 


0.025 


150 


E 


3 • 10" 4 


IOO/55 


0.5 


1.045 


0.31 ±0.02 


21 


32 


0.025 


150 


Q 


2 • 10" 2 


IOO/55 


0.5 


1.0648 


0.409 ± 0.003 


22 


12 


0.1 


2.614 


E 


9 • 10" 8 


/22 


1 


1.0052 




23 


32 


0.1 


10 


E 


4 • 10~ 6 


10/ Bo 


1 


1.0260 




24 


32 


0.1 


10 


E 


8 • 10~ 8 


He 


1 


0.7266 




25 


32 


0.1 


10 


E 


4 • 10~ 5 


10/5a + fa 


1 


0.9694 





Table 1: Runs performed. Here At represents the integration step, T the total integration time and e the relative 
numerical error on the constrained quantity Q m . The following notation is also used: fa = fa + fa, with 
fa = 0.00015 on k = (4, -3) and fa = 0.000H on k = (3,4); fa = f 2a + fa, with fa = 0.0006 on k = (2,0) 
and fa = 0.0004i on k = (0, -2); fa = -0.00009 + 0.00008i on k = (0, -2). Values of {Pi, m )/vi which agree with 
the (extended) EC within 15% error are marked in boldface, where it is intended that an error may affect the last 
digit. The values c 4,m and the associated errors were computed by fitting the final part of curves like those in Figj6| 
with horizontal lines, when sufficient statistics was available. 
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Figure 2: Run #11: behaviour of E, Q, P and f3\^ in the NS "buildup" part and once P is constrained. The 
format of presentation is the same of Fig||. Also here is close to v\. 



quite remarkable if we consider that the only information we pass from the NS to the GNS system is the NS field 
at an arbitrary time t after the stationary state has been reached, and that two consecutive snapshots may yield 
rather different values of the chosen Q m - In doing so, we assumed that any phase space point of the steady state 
trajectory carries enough information for the GNS run, and our results validate this approach. 

The fact that good agreement is obtained when imposing the conservation of a quantity like E or Q is not 
entirely surprising, and might be connected with the notion of "rugged invariants" , i.e. of global quantities which 
are approximately conserved in spontaneous turbulent decay. While the constancy of such quantities has been 
imposed for structural reasons here, that was earlier the key assumption of a number of theories of turbulence. 
Without attempting to make a review, we remark that the stationarity of E, together with that or not of Q, were 
the basis of various models of either stationarily forced or freely decaying two-dimensional turbulence. Among 
these we quote: a) the equilibrium statisical mechanics of Fourier modes, with fixed E and Q, due to Kraichnan 
|p5| ; b) the minimum enstrophy theory, in which the decay is ruled by a minimization process of Q at fixed E |p6f ; 
c) the maximal entropy theory, in which the decay involves again a lowering of Q at constant E, but the decrease 



in Q is explained probabilistically as loss of information due to mixing 27 ■ In our cases, however, the hierarchy 



of conserved quantities implied in these theories seems to be absent: the equivalence between the stationary state 
of the GNS and of the NS equations is found for all constrained Q m . On the other hand, the amplitude of the 
fluctuations of /3/ jln increases with m, as can be seen comparing Figs.|l| and | 

Another interesting fact is that in the NS equations we can fix the value of / and of vi independently of each 
other, while the amplitude of the fluctuations of [ii^ m grows as the ratio f/ui increases (see Fig.||). This is different 
from the case of particle systems, where the dissipative term is proportional to the forcing term. Also, differently 
from the case of particle systems, the quantity uik in the cut-off GNS equations has k 21 Pi^mUJk instead of just (3^ m uJk- 
It is striking, then, that the equality between (f3i >m ) and v\ is verified irrespective of the value of /, as seen varying 
/ over several orders of magnitude. 

To check the validity of the fluctuation relation Eq. (p7|) , we considered a number of long runs, which we cut in 
segments of given length r, for various values of r, for different realizations of the forcing term, and for different 
values of the viscosity v\. Moreover, we considered four different constraints in the cut-off GNS equations: E, Q, P 
or H ^constant. Sample histograms of the distribution of the averages of (3 1 _ 1 over the times r are shown in Fig]4] 
for constant E and a few different values of r. Clearly, the distributions become more and more peaked around the 
mean (/3i,_i) as r is increased (large fluctuations of /3i,_i are wiped out by the average). Such distributions deviate 
however from gaussian, and the non-gaussian character is not mildened by the averaging; the actual computation 
shows that the kurtosis ((/3; jrn — l) 4 )/((A jTO — l) 2 ) 2 is significantly different from 3 (the value for the gaussian 
distribution), and that this difference seems to increase both with r and / (see Fig||). This is not exactly the 
case for particle systems (at moderate forcings), whose probability distributions can be interpolated quite well 
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© M=32. F=f5S (nm#l) -7 M=32, F=10O*f55 (run#9) 
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Figure 3: Histograms of /3i,-i with error bars, for runs with constrained E, original v\ — 10 4 and various forcings. 
Note the coalescence of histograms referring to runs at different resolution but equal forcing. 



by Gaussians, although they should not be Gaussian distributions |13[. Similarly to previous works [£8|, |l5| for 
particle systems, we do not separate in time the adjacent segments to decorrelate them. The reason is that, in all the 
cases we considered, the only effect of the decorrelation was the worsening of data statistics. Our test consisted in 
computing the left hand side of Eq. ( p7[ ) for many values of p, in verifying the consistency of the results with a linear 
law of slope C^' m (r), and in extrapolating C^ m (r) to large r. Given the observed proportionality between m and 
tf£ s , Eq.@ holds for the distributions of a^ s , with the slope C^ m (r) replaced by C l f n {r) = C l f > (t)u 1 / 'a? TS '. 

Our results show that the distributions of the values /3 ; m are in all cases consistent with the linear law Eq.(p7|). 
The fit to a straight line, performed according to standard least-squares minimization, is always remarkably good, 
as in the cases shown in Fig.^. The error bars are derived from the errors on the histograms, which, according to 
standard statistical analysis, are estimated as the square root of the bin count. The value of the linear regression 
coefficient vs. r is plotted in Fig.|| for several cases. The slope of these lines decreases with t, and converges 
to a finite value at large r when E, Q, P or H are fixed. Therefore, the fluctuation relation is valid in these 
cases, although presently we cannot compare our results with Gallavotti's theoretical prediction || which links 
c l f n = lmv^oo C p m (t)vi / s to the Lyapunov spectrum of the system at hand. However, consistently with these 

predictions, we obtain < c^ m < 1. We conclude this section noting that the dependence of C^' m (r) on r (cf. 
FigJ^) is different from that observed in the work fll|| , while it is similar to that observed in the particle systems 
ofRef.lil. 



5 Conclusions 

1. We investigated the equivalence between the NS and GNS equations. Such equivalence is verified for a variety 
of cases, and actually appears more robust than expected in B. Indeed, (/3i, m ) approaches v\ at least when E or 
Q or P or H are fixed in the cut-off GNS equations (cf. Table p. Moreover, neither the thermodynamic limit nor a 
large number of modes is needed to obtain this result. This leads us to conclude that the EC can be extended to the 
case of constrained Q m , for any m > — 1. The equivalence of the NS and GNS equations can be equally expressed 
in terms of (J3i,m) an d y h ° r 01 (°f ? m S ) an d cr / VS ! due to the proportionality between of^ 5 and /3z jm expressed by 
Eq.(|23|). We preferred to use (/3i. m ), which is directly connected with a quantity of ordinary hydrodynamic interest: 
the viscosity of the fluid V\. Moreover, the values of c^ m appear to cluster around a single value, independent of 
the resolution of the run and of the parent viscosity, as can be inferred from Table [j] using c^ m = c vija^ s . This 
is not the case for c^ m . 

2. The fluctuation relation Eq.(p7|) is seen to hold when E,Q,P,H are constrained. In these cases, the left 
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Figure 4: Normalized histograms of l for various values of r, for the run #9. The values of the kurtosis evidence 
the non-Gaussian nature of these distributions, which becomes more and more pronounced as r grows. 
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Figure 5: Linear fit of 
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to (7g 1 (r) p for two different averaging times r, from run #9. 




Figure 6: Convergence of the slope C,J' m (r) of the regression line for various runs. The legend stresses the resolution 
M, the quantity constrained, the force used and the parent viscosity, also reported in Table [|. The data presented 
here refer respectively to the runs #11, #10, #20, #5, #12 and #14 of that table. Note the overlap of the curves 
of runs #10 and #20. 



hand side of Eq.(|7|) is perfectly fitted by straight lines at all values of r, (see e.g. Fig.||) with slopes converging 
to definite values c^ m in the limit of large r (cf. Fig^j and Table Q). This at one time validates the CH in the 
framework of hydrodynamics, and yields a new result related to the NS equations. Moreover, the values c^ m are 
consistent with Gallavotti's predictions (cf. Eq.(4.1) of g), which imply < c l ^ m < 1. 

3. Our results concern statistically stationary states of the dynamics which, from a fluid-dynamicist point 
of view, are achieved at very late times. The resulting flows are substantially dominated, in absence of a proper 
large-scale dissipation mechanism, by vorticity structures of the size of the computational domain. These large 
scale structures are seen to be unsteady, fact which is by itself responsible for the fluctuations of the system. In 
such conditions the mechanism of transfer of energy toward small k has already pushed almost all the available 
energy to the largest accessible scales; the formation of a self similar direct enstrophy cascade or of an inverse 
energy cascade are prevented. If one is concerned with universal inertial range properties of the turbulence, such 
kind of states would be inappropriate because they are dependent on the finite size boundary conditions, hence 
not universal. Our point of view is indeed different, and we do not consider the poor resolution or the lack of a 
universal turbulent cascade as shortcomings. We are interested in the statistical behavior of a given dynamical 
system, and the computational limitations on N do not hinder our investigation. In this spirit, we treated the 
viscosity coefficient as a mere parameter of the simulation; we did not match resolution and viscosity in order to 
guarantee the resolution of a dissipative range. Again, the point is that we investigate special properties of the 
vorticity modes which, because of their dynamical nature, should not depend on the level of resolution. 

4. Some confusion is sometime made on the nature of the GCFT of Ref.Q, because the p.d.f. obtained for 
the fluctuations of the entropy production rate in particle systems are indistinguishable to the eye from Gaussian 
distributions. This led some to believe that the statement of a linear law for the odd part of the p.d.f. of the 
fluctuations in particle and hydrodynamic systems is a quite generic result and just a manifestation of the validity 
of the central limit theorem. This belief, however, is not well founded. In the first place, the relation investigated 
here and the GCFT concern large and not small deviations, hence Gaussian distributions should not be expected. 
For instance, these distributions have finite support for m = — 1, 0, because the values of f3i iTn are bounded, being a 
ratio between linear and quadratic functions of the vorticity field. Furthermore, differently from ]l3| , the kurtosis 
of our distributions are manifestly different from the Gaussian value 3 (cf. Fig.||) at finite t's, and we haven't 
observed any convergence to 3 for growing r. Even in the case of in which a Gaussian interpolation of the 
data is possible, a connection with the central limit theorem is far from obvious. 



5. Two recent papers, [[[T| and considered the verification of the EC and of the fluctuation relation in fluid 
dynamics, suggested by Gallavotti [pL The present work strengthens these earlier findings, along with new results. 
In fTT) , one experimental fluctuating time series, representative of a heat flux, is analyzed. The p.d.f. of this series 
is evidently non-gaussian, but its odd part can be fitted by a linear relation, as in our cases. The asymptotic value 
of the slope of the fitting lines is not given. In contrast, we were able to explore several cases and to check how the 
asymptotic c^ m are reached. 

The work of , instead, is mainly concerned with the relation between ordinary shell models and reversibly- 
damped counterparts. The authors find differences in the statistics produced by the two versions of the shell 
models, as well as limitations on the values which the parameters can take for one kind of equivalence to hold. 
Such equivalence is found in the correct reproduction of the energy cascade and of the multiscaling of the structure 
functions. No analysis of the GCFT is attempted there, at variance with our paper. It is however interesting to 
compare our and their approach, as it seems that the dynamical equivalence between "normal" and "reversibly 
constrained" dynamical systems representing fluid turbulence is not always as complete as in the NS-GNS cases. 

6. One final remark concerns the debate on the connection between the GCFT and the so called Evans-Searles 
identity (ESI) Jl(| |2jJ. The ESI is a relation which concerns time reversible dynamical systems, and the Liouville 
mesaure on the phase space of such systems. Hence it can also be applied to our truncated GNS equations. 
In particular, let E p C f2 be the subset of initial conditions of trajectories along which the phase space contraction 
after an evolution of any length T is e^ va< ^- m T . Then, the ESI implies that |^9| 

log n L (E p ) - log fi L (E. p ) = p(<t^ s )T . (28) 



This equation is formally similar to Eq.(27) re- written in terms of <jf^ s , for T = r and C^. ,m (r) = 1. Now, the 



fluctuation relation of the GCFT also has a slope of 1, in the restricted domain of dynamical systems with dense 
attractors not considered here. This led some to believe that the ESI and the GCFT describe the same quantities, 
namely the fluctuations of the phase space contraction. However, our results explicitly show that the ESI cannot 
describe these fluctuations, because we have c 1 ^ 171 < 1. The ESI, instead, concerns the relative probability of 
independent "trajectory histories" emanating from the Liouville distribution. 
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